NA5A IN D-861 


WFf* % <?*"'** rt I . * ^ ;*■ % ?» * , Ta 

p*^ I ir'' ^ ^ ^ ^ ^ 



C / 

NASA TN D-861 


LIBRARY COPY 

SEP 20 1961 

SPACE FLIGHT 
LANGLEY FIELD, VIRGINIA 


TECHNICAL NOTE 

D-861 


f 


I 


AN APPROXIMATE ANALYSIS OF FILM COOLING ON BLUNT BODIES 

BY GAS INJECTION NEAR THE STAGNATION POINT 

By Byron L. Swenson 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON September 1961 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D- 86 l 

AN APPROXIMATE ANALYSIS OF FILM COOLING ON BLUNT BODIES 
BY GAS INJECTION NEAR THE STAGNATION POINT 
By Byron L. Swenson 

SUMMARY 


An approximate method for the estimation of laminar heat transfer 
to ‘blunt bodies with gaseous film cooling is developed. Attention is 
focused on the parameters which are important for the design of an 
attractive heat protection system. Application of the analysis is made 
to calculate the approximate coolant weight requirement for both a 
circular and a parabolic entry. 


INTRODUCTION 


Continued study of heat-protection systems has produced a number of 
efficient means for dissipating or absorbing the aerodynamic heat other- 
wise convected to a re-entry vehicle. Among the systems which have 
shown promise are those based on the injection of foreign gases into the 
boundary layer of a vehicle. In particular, it is well known and under- 
stood that transpiration or normal injection of a foreign gas from a 
porous wall can efficiently reduce the local heat-transfer rate (see 
ref. l). This transpiration method does have some disadvantages, however. 
First, the injected gas reduces heat transfer primarily in the immediate 
area of injection and has little influence on heat transfer downstream 
of this area (ref. 2 ). In addition, a porous wall presents some structural 
problems. Since these disadvantages are associated primarily with the 
method of injection, it seems proper to inquire if other means of injec- 
tion can be used which are more attractive structurally and which hold 
promise of providing a larger over-all cooling effect. One such method 
appears to be the use of gas film cooling. With this method the gas is 
injected tangentially to the wall and usually from a single port. With 
tangential injection, mixing of the injected gas is reduced and the gas 
tends to be maintained in a thin film. 

Film cooling has been used to advantage in the heat protection of 
both wind-tunnel nozzle walls (ref. 3) and blunt re-entry shapes (refs. 

4, 5, and 6 ) . In addition, it has been shown that the streamwise gas 
injection scheme can be used to reduce the skin-friction drag of slender 
bodies (ref. 7) • 
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For the most part; the previous investigations of film cooling have 
been experimental in character. While these investigations have indicated 
attractive features of film cooling; it would seem that the virtues of 
the method might be further exploited if a simplified analytical treat- 
ment of the phenomena were available. These considerations then define 
the primary objectives of the present paper. The first objective is to 
obtain an approximate analytic description of film cooling and the second 
objective is to employ the results of this analysis to determine the 
design criteria which make the method most attractive. 


SYMBOLS 


A reference area, ft 2 

B distribution parameter 

b span width, ft 

Cp drag coefficient 

Cf skin- friction coefficient 

Ci injection coefficient 

Cp specific heat at constant pressure, Btu/slug °R 

Cp pressure coefficient 

AH enthalpy potential, cp(T r - %) , Btu/slug 
h heat-transfer coefficient, Btu/ft 2 -°R 

Ip drag surface integral 

Ig weighted surface area integral 

k thermal conductivity, Btu/ft-sec-°R 

m-p mass rate of film flow, slugs/ sec 

m-p total mass of film gas required for protection during entry, slugs 
M molecular weight 

P body perimeter, ft 

p pressure, lb/ft 2 
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q_ heating rate, Btu/ft 2 -sec 

Q integrated heating rate, Btu/sec 

r "body radial coordinate , ft 

R characteristic dimension of body (nose radius) , ft 
R universal gas constant 

s arc length along body surface measured from the injection port, ft 
T temperature, °R 

t time, sec 

ti time at which gas film cooling is turned on, sec 

tz time at which gas film cooling is turned off, sec 

u velocity, ft/ sec 

V vehicle flight velocity, ft/ sec 

W vehicle weight 

rj_ injection coefficient 

& local body slope, radians 

5 C cone half -angle, deg 

by. wedge half -angle, deg 

p. coefficient of viscosity, slugs/ft-sec 

p density, slugs/ft 3 

Subscripts 

0 condition of film at point of injection 

1 air-film interface 

A air 

b evaluated at the base 


F film 

H helium 

r recovery condition 

s stagnation value 

w wall condition 

m ax maximum value encountered during entry 
oo free-stream condition 


ANALYSIS 


As noted earlier, the following analysis is designed to provide a 
simple description of film cooling. In keeping with this objective, 
only laminar flows will "be considered and it will "be assumed that "both 
the film and air are ideal gases with Prandtl numbers of unity. From 
the analysis point of view, the problem can be treated in two parts: 
first, the case where the film is very thin compared to the air boundary 
layer and, second, the case where the gas film replaces at least the 
boundary- layer thickness of the air flow. The analysis of the two cases 
is as follows. 


Very Thin Film 



Local heating rate .- The 
thin- film model for film cooling is 
predicated on the assumption that the 
gas film is very thin compared to the 
m F air boundary layer. This abstraction 
of the actual boundary- layer flow is 
shown in sketch (a) where typical 
velocity and temperature profiles are 
shown schematically. For simplicity 
it is assumed that there is no mixing 
across the air-film interface , or, by 
implication, that Ui is the same as 
existed in the air boundary layer a 
distance bp from the wall without 
injection. It is also assumed that 
Ux < < u 


Sketch (a) 
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With the very thin film assumption in mind, the heating rate to the 
wall can be expressed as follows: 



Ti ~ T v 

Sp 


( 1 ) 


The thickness of the film can now be established from considerations of 
the mass flow in the film and the wall shear stress. They are, 
respectively: 


mp = CpiUi&jpP 


( 2 ) 


where 


P 


f2jtr( s) 


axi symmetric flow 
planar flow 


C 



&F “*• 0 


and 

= \ PA u A 2c f A (3) 

Attention is directed 
for the moment to the vis- 
cosity in equation ( 3 ) . 

For simplicity it will he 
assumed that the viscosity 
may he approximated hy 
p = kT where k is the 
same for all gases. (The 
accuracy of this approxi- 
mation is indicated in 
sketch (h).) This approx- 
imation is perhaps hest for 
gases , such as air, helium, 
nitrogen, and oxygen, and 
it is poor primarily only 
for hydrogen. In any event, 
the approximation is adopted 
here. Assuming no pressure 
gradient normal to the wall 
through the boundary layer 
we may then write the 
following: 



Sketch (b) 



where M is the molecular weight. 


Equations (2) and (3) together with equation (4) can now he solved 
for the film thickness; thus, 

m a . V/2 

Mj? m F 


&F = 


c fA? 




Ta \ PA U A 


Substituting equation (5) into equation (l) using 

kp 1 


hf^i c : 


Vj? 


= 1 


(5) 


yields 


VaP 


% = \ PA u A c f fi^ 1 ! 


% 


1/2 


(6) 


where the enthalpy potential is AHi = c Pf (T! - T w ) . By Reynolds analogy 
the heating rate without a film is given hy: 


!A = \ PA u A c f A ^A 


therefore, 


/ Mp \ ^ 2 
q p / PA P ^ \ ^ 

1 A l *F c f A ) ^A 

It can he seen from equation (7) that a film of a low molecular weight 
gas with a low enthalpy potential would reduce the heating rate. 

Heating-rate distribution .- The distribution of the cooling effect 
is determined hy a heat balance on the film. With reference again to 
sketch (a) , it is assumed that the amount of heat conducted across the 
interface is approximately that convected to the wall in the absence of 
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the film. This assumption appears reasonable if Ti < < T ma x . With the 
additional assumption of an isothermal wall, the heat balance is expressed 
by: 

(q_A - q.p)P ds = — m F dAHi 

The result can be rewritten into the following form with the assumption 
that for blunt bodies in hypersonic flow AHa is constant over the body; 




lit 


, AHn 
a 


W <lA,s R R 2< lA,s R2 ^A 


(8) 


where R is a characteristic length of the body (say the nose radius). 
The heating-rate distribution is now determined as a function of body 
shape, mass rate of injection, and flight conditions by solution of the 
linear differential equation resulting from the combination of equations 
(7) and (8) . 

+ 2 9 a , s^ p %A P ^1 _ 2 1 a,s r2 ^A P 
&(s/R) mpAH^ ^A,s R ^A mpAHA ^A,s R 


The solution is as follows: 

< 1 F _ 3 


where 


«A E( s/R) l(^1) +F (r 


P - 


p A p(MjyM A ) - 

l/2 

r(Mp/M A )AH A [i A P A u A Pl 

. & F c f A j 


L 2^,b 9a/1A,sJ 

a.sR 2 

; ps/R 

n lA P Jf ,\ 

exp — ; L 

\ mpAH A 

Vo 



- 1/2 


( 9 ) 


A.iAsi r /E 


R / 


E(S) 

<U,s R 


and AHx/aHa is the enthalpy-potential ratio at the injection point. 
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The parameter (3 can he given physical significance if the equation 
for the film thickness (i.e., eq. ( 5 )) is rewritten in terms of (3. 

5 = u> Vi 

F A Ta (l/2)p A u A 2 c fA P 


or 



Comparison of this result with equation ( 3 ) shows that (3 can he 
interpreted as 


(3 ~ ua/ui 


Since it is necessary that ua/ux "be much greater than 1 to satisfy the 
approximations of the analysis, it follows that (3 > > 1 is a requisite 
of the analysis. It also follows that a restriction must he imposed, on 
the mass flow in the film so that (3 > > 1 is satisfied; that is, 


mp < < 


p a p(m f /m a ) 

C *A 


(m f /m a )ah a P A P A u A P 
2 ^A,s V^s 


( 10 ) 


The results obtained from equation ( 9 ) for film mass flows in excess of 
the above restriction will he questionable since this condition will 
violate the basic equations of the analysis - in particular A equation ( 3 ) . 

Integrated heating rate ■ - The heating rate integrated over the body 
can not be obtained in closed form for this type of film flow. However , 
it may be obtained by graphical integration of equation ( 9 ) • 


Thick Film 



Local heating rate .- The basic 
assumption of the thick-film model of 
film cooling is that the gas film replaces 
at least the boundary -layer flow of the 
air. Again mixing across the air-film 
interface will be assumed small A being 
only on a molecular scale. This abstrac- 
tion of the flow is shown schematically 
in sketch (c) . It is obvious that there 


Sketch (c) 
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can "be no pressure difference across the interface and that the assumption 
of smal 1 mixing requires that at the interface, the velocities of the air 
and the gas film be equal. 

With these points in mind, the ratio of the local heating rate with 
film cooling to that without can be written as 

Ip hp (T r - T w )p 

0 .A -^A ^r ” ^w^A 



Also, by Reynolds analogy, the ratio of the heat-transfer coefficients 
is 


h F (pUCfC p ) F 
hA “ (puc f Cp) A 


Combining equations (ll) and (l 2 ) yields the following expression: 


q F (pucfOjAHp 

qj “ (puc f ) A AH A 


where 


AH — Cp(T r - T^) 


The densities in equation (13) are evaluated with the equation of 

state 


P = 



T 


where R is the universal gas constant. With this result and with the 
previously noted conditions, that is, up/u A = 1 and pf/pa = 1> equation 
(13) becomes 

% ^A / , n 

q A _ M A Tp c fA AH a 

Attention is now directed to the ratio of skin-friction coefficients. 
This ratio was evaluated in reference 7 for the special case of a helium 
film. The result of reference 7 can be generalized to other gases if 
the use is again made of the approximation that viscosity may be written 
p = kT where k is the same for all gases. Finally, since the 


I 
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skin-friction coefficient is inversely proportional to the square root 
of the Reynolds number, the following ratio is obtained: 


C ?F _ ^F /% 

c f A " t A n/% 

Substitution of equation ( 15 ) into equation (l4) yields: 

q_F I^F 
OA ” J MA £§a 


( 15 ) 


(16) 


From equation (l6) , it is again seen, as in the previous flow model, 
that low heating rates will he obtained with a film of gas having a low 
molecular weight and a low enthalpy potential. 

Heating-rate di stribution . - With the ratio of local heating rates 
determined, consideration must now be given to the distribution of heating 
rates over a body with film cooling. To obtain this distribution, it is 
necessary to evaluate the enthalpy potential of the film, AHy. The 
potential increases as the film flows aft over the body as a result of 
the absorption of heat by the film. As the enthalpy potential increases 
so does the ability of the film to convect heat to the wall. 

A limiting estimate of the distribution of the enthalpy potential 
of the film is immediately evident. The simplest assumption is that the 
enthalpy potential of the film does not change as it flows around the 
body. This would occur if the film not only replaced the boundary layer 
but was sufficiently thick that conduction and diffusion of heat from the 
exterior air did not penetrate to the boundary layer by the time the film 
flowed off the rear of the body. In this case, equation (l6) would 
completely describe the heating -rate distribution about the body since 
AHp would be constant at the value for injection. Unfortunately this 
condition will not often exist in practice, especially for axi symmetrical 
bodies where the film generally will thin greatly as it flows aft. 

Integrated heating rate .- In the case of a thick film it is easy to 
determine the integrated heating rate or integral of the heat rate over 
the body surface. The heat accepted per unit time is 


Q = 


% 

qP ds 


where is the arc length to the base of the body. It is easily 

shown that the ratio of the integrated heat rate with film cooling to 
that without can be written as 
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F^F.o 
Q a 7 Mft AH a 


( 17 ) 


The term AHp 0 is the enthalpy potential of the film at the injection 
point. 


Intermediate Film Model 


A third possible model of the film flow has characteristics somewhat 
intermediate to those of the two models already considered. With this 
model; it is assumed that the film replaces just the boundary- layer 
thickness of the air flow. 

Local heating rate .- With this model; the local heating rate is 
again given by equation (l6) . The difference between this model and the 
previous one lies in the amount of heat transported across the air-film 
interface and; hence; in the distribution of heating rate. 

Heating-rate distribution .- A limiting condition on the amount of 
heat transported (i.e. / by diffusion and conduction) to the film from 
the exterior air across the interface is that it is the same as that 
convected to the body in the absence of a film. At first glance; this 
assumption appears to be contradictory to that used to develop equation 
(l6) . However; it is felt that; especially when a large potential exists 
between the film and free- stream conditions; the conduction across the 
interface could be of the same order as the convective heat transfer. It 
should be noted that a similar assumption was employed, by the authors of 
reference 8 to obtain a correlation of recovery temperature for film- 
cooled walls. This assumption appears to be conservative in the sense 
that less heat will probably cross the interface; but this does not mean 
that the entire analysis is conservative. There are effects not considered 
here; such as mixing and diffusion; that may have a detrimental effect 
on the cooling provided by the film. In any event; due to the simplicity 
gained from this assumption; it will be adopted here. 

A heat balance of the film based on these considerations is given 
by the following. 


(q.A " Ip) P ds — m.p d(AHp) (l*8) 

Equation (l8) can be rewritten in the following form with the assumption 
that AH^ is constant over the body. 

/ <lF\ <1 A p s “F^A ^F 

V “ W <3A,s R R ~ ah a 


( 19 ) 
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The heating-rate distribution with film injection can now he determined 
as a function of body shape, mass rate of injection, and flight conditions 
by the solution of the linear differential equation resulting from the 
combination of equations (l6) and ( 19 ) • 



+ 



mpbH/v 


JA_ 

9A,s 


p AHp 


1A , s r2 P 

mpAH A q^g R 


The solution is as follows: 


IF , ( ^ /"% ^FjO^ -B 

q A = 1 A W MA AH A ) 6 


where 



and 



r s/R <3a P s 

Jo As R R 


(20) 


The term AHp Q is the enthalpy potential of the film at the injection 
point . Since 


1 


fi F %,o _ / 

J M A - AH a V 1 - q A / o 


equation (20) can be placed in a more symmetrical form 

1 - (if/ia) _b 
[1 - (qp/lA) ] Q 


( 21 ) 


In the discussion that follows, the term [l - (lp/l A )l 0 will be called 
the injection parameter and the term B, the distribution parameter. In 
addition, the integral I s may he thought of as surface-area integral, 
weighted as to heating rate. 
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Integrated heating rate .- In the case of the intermediate-film model 
it is again easy to determine the integrated heating rate or integral of 
the heat rate over the "body surface. This integral helps define the 
total heat load accepted by the body. The heat accepted per unit time 
is again 


, s b 


Q = 


qP ds 


Therefore, the ratio of the integrated heat rate with film cooling to 
that without can be written: 


Qp 

q a 




1A <lA,s R R 

' S ^ R P s 

3 *A,S R R 


or 


Qp 

Q a 


plsCsb/^) Q.-C, 

I JL (I s ) dl s 

o 


I s (sb/R) 

Sub st itut ing equation (20) into equation (22) and integrating yields 


(22) 


Qp 

1 " Qa 


i - 


q.p\ i _ 


-Bb 




( 23 ) 


where is evaluated at the base and is 


*b = 


'% pS-p fb 

Ma q ‘A,s iri s R 


ApAH A 
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RESULTS AND DISCUSSION 
Comparisons With Existing Experimental Data 


In an attempt to assess the accuracy of the present analyses, 
comparisons will he made with existing experimental data for the effects 
of film injection on heat transfer • Some of the available data on film 
cooling are presented in reference 4 where results were obtained with 
both helium and nitrogen as the injected gas. For these gases results 
are presented in the form of the integrated heating rate to a hemispherical 
nose for segments extending 30°, 60°, and 80° from the stagnation point. 
These experimental results are presented as a function of the injection 
coefficient, C-j_, in figures 1 and 2 . 

Additional data on film cooling for the case of air injection are 
obtainable from reference 6. Again, only integrated heating rates to a 
full hemisphere were obtained. For the purposes of comparison with the 
results of the present analysis, the data of reference 6 have been con- 
verted from the mean heat-transfer coefficients presented to the form of 
integrated heating rates. This conversion has been made first with the 
use of the measured recovery and wall temperatures and secondly with the 
use of the measured recovery temperatures and assuming a cold wall. The 
latter conversion, of course, assumes that the measured heat-transfer 
coefficients are independent of wall temperature . It should be noted 
that the approximations made in the analysis will be most closely satis- 
fied when a large potential exists between the wall and the free-stream 
recovery conditions . Since the data of reference 6 were obtained with 
the use of a heated model with a small temperature potential between 
the wall and free-stream recovery conditions, it seemed necessary to 
adjust the data of cold wall conditions before comparing them with the 
results of the analysis. The data from reference 6 are shown in figure 3. 

For each set of data, comparisons are made with the results of each 
of the present analyses. The theoretical estimates shown in figures 1, 

2, and 3 were obtained with equations (9), (17), and (23) , In each 
application of the thin-film analysis it was found that the restriction 
placed on this analysis, as given by equation (10) , was seriously violated. 
Typical values of p were about 1 / 4 , indicating that the velocity at 
the film-air interface was four times the local free-stream velocity. 

It is suggested, therefore, that for the available experimental results, 
the mass flows were too great to permit valid use of this analysis. For 
this reason, the reasonable agreement obtained in some cases between the 
results of this particular analysis and the experimental data must be 
viewed with some reservations . 

Perhaps the most obvious conclusion from the comparisons shown in 
figures 1, 2, and 3 is that the thick-film analysis is the least accurate. 
Apparently this analysis does not represent the film-cooling phenomena 
except as a limit as the mass flow becomes very large. 


15 


The results obtained with the intermediate model appear to represent 
a reasonable approximation for most of the experimental data. Detailed 
examination of the three sets of data indicates some inconsistencies , 
however. The helium- injection data of reference 4 seem to be quantita- 
tively predicted, at least at the lower injection rates, by equation (23) . 
The trend of the data of reference 6, when adjusted to cold-wall condi- 
tions, is also predicted by this equation. On the other hand, the 
nitrogen-injection data of reference 4 disagree with equation (23) of 
the analysis. Essentially no reduction in heating rate was measured in 
reference 4 at the lower injection rates. It should be noted that this 
result is not in agreement with the results of reference 6 for a gas of 
approximately the same molecular weight . 

In any event, it would appear that the available data do not permit 
any definite conclusion as to the accuracy of the present analysis. 

Real justification of the analysis and the conclusions gained from it 
will have to wait for additional and more comprehensive film- cooling 
data. In spite of this restriction, further considerations of the 
application of film cooling are warranted. For these further considera- 
tions, attention will be restricted to the intermediate-film model. 

This model will be used since it is at least as accurate as the others; 
it is not subject to the reservations attached to the thin-film model; 
and it is somewhat more simple in form. 


Applications of the Analysis 


Injection coefficient .- With the dependent effects of film injection 
defined in the analysis presented earlier, attention at this point is 
turned to the independent variable, the rate of injection. In most 
studies, for the purposes of scaling, the rate of injection has been 
described in terms of an injection coefficient, C^, which may be written 
as 


r. = mF 

i p u jtR 2 


This coefficient refers the injection rate to the rate at which free- stream 
air is swept out by an area of radius , R. While this coefficient is 
convenient, it holds no true physical significance. The results of the 
analysis of the intermediate -flow model suggest the use of another 
coefficient of injection which appears as a component of the parameter, 

B: 


m-p 

r i = (^TSftjP 




( 24 ) 
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When this coefficient is used, the rate of injection, m F , is referenced 
to mass rate of coolant it would take to at sort the heat convected to a 
stagnation area of R 2 if that mass were increased in enthalpy ty an 
amount equal to AH/^. 1 For this reason, the coefficient would seem 

to have greater significance for correlations than does Cq. Thus the 
distritution parameter, B, may te written as 


B = 



(25) 


Coolant requirements .- The total weight of film coolant used during 
a given flight may he written 


nip = / m F dt 
^ t 

where this equation is integrated over the time of flight. In terms of 
the injection coefficient, Fq, this equation is 

or in terms of the distribution parameter, B, it is 


mp = 


It is assumed in equation (26b) that for hypersonic flight conditions, 
the ratio q A /q A s is independent of Mach number and thus I g is not 
a function of time . 

Film injection may be applied to provide cooling in several different 
ways. For example, it may be desirable to inject in such a way as to 
keep the integral of the heating rate over the body at some constant 
fraction of that with no injection during the entire flight. For this 
case and for AHp^o < < AH^ then Tq or B is a constant (see eq. (23)) 
and not a function of time. Therefore, 

1 Wote that the reference mass rate differs by a factor from the mass 

rate of coolant it would take to absorb all the heat normally transmitted 
to the body if that mass were increased in enthalpy by AH/^. The factor 
is I s (s t /R). 
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mp = 



q A ; s 


dt 


In some cases it might also he desirable to inject in such a way 
as to keep the local heating rate at any point on the body below a 
certain design limit. In this case then, equations ( 2 l) and (25) must 
be used together with the known heating-rate distribution of the body 
without injection to determine how should vary with time along the 

flight trajectory. The total mass can then be calculated by the use of 
either equation (26a) or (26b) . 


Another possible application of film cooling might be to keep the 
integrated heat rate below a certain design limit. Consider, for 
exanple, an L/D = 0 entry with an integrated heat rate as a function 
of time as shown in sketch (d) . Consider now a structure with a certain 
design integrated heating rate, Qp, that it can withstand or dissipate 
Without assistance. Obviously optimum performance of the heat protection 
system is attained ideally when 
the structure is used to its 
fullest capacity. Therefore, 
the film-cooling system would 
not be turned on until time t^. 
and then would be controlled in 
such a way that Qp = Qp until 
time t 2 when the heating rate 
Qa would again become equal to 

Qp. 

To keep Qp equal to Qp 
during the time interval tj. to 
t 2 it is necessary to control 
the mass rate of injection (and 
hence B) as indicated below. 

From equation (23) 



Sketch (d) 


Qf _ Qd / _ <2.f\ 1 - e“ Bb 

Qa Qa V W 0 B b 


Consider now injecting the film so that £Hp^ 0 = 0 ; then the parameter 


1 _ IF 


1 


o 
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thus 


Bb(t) 


Q A (t) 


1 - e "®b(t) “ Qp 


(27) 


A graph of the left side of equation (27) is shown in sketch (e) . Suppose 

that for example Qp = (i/2)QA,max, 
then must vary over the period 
ti to t s as shown hy the arrows in 
sketch (e) • It is easy to see that 
during this variation of 





- e-Sb 


B =b 


Thus we may write (noting that 

= ^a) j 


B b 

Sketch (e) 






. ft) . ajh ' A '~' - 

^ %(t)AH A (t) Q A (t) - Qp 


or that 


“p(t) s . 


1 Q A (t) - Qp 


% AH A (t) 

Equation ( 29 ) can he now rewritten in the form 


( 28 ) 


(29) 



V\ q-A,s (t) 
.V AH A (t) 


1 


% " I 

Q A (t)_ 


or the total coolant mass required is 



( 30 ) 


(31) 


These are not the only ways in which the coolant injection might 
he controlled during given flight conditions. It is usually easy to 
determine how the distribution parameter and hence the mass rate 
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of injection, must vary with time for the particular application under 
consideration. Then either equation (26a) or equation (26b) can be used 
to calculate the total mass of coolant required. 


Design Parameters 


Let us now inquire what information can be obtained from the 
intermediate-film model as to the nature of the film and body shape 
which will provide the heat-protection system with the best performance. 
Attention will be given first to the most attractive type of film. 

Injection parameter .- The present analysis indicates that the 
reduction in local and integrated heating rates is directly proportional 
to the injection parameter 

/ 9p\ _ S m F,o 

V " Vo “ n/ u a a% 

It can be seen from this result that the cooling effect of the 
film increases with decreasing molecular weight, Mp, and with decreasing 
enthalpy potential, AHp t0 , of the injected gas. Thus a light-weight gas 
injected at low temperature would appear to be most attractive for use 
in film cooling. This result is, of course, in agreement with the con- 
clusions reached in previous studies of the effects of film injection 
(see, e.g., refs. 3> and 7) • The body shapes which appear 

attractive for film cooling are, in part, dictated by consideration of 
the cooling distribution. 

Distribution parameter .- The distribution of the local and integrated 
heating rates provided by film injection depends upon the distribution 
parameter, B, where 


_ n/ Mp/M^ qA.,s R2l s 

mjAH A 

This dependence is defined by equations (2l) and (23) and is shown 
graphically in figure 4. It is apparent from figure 4 that from the 
heating standpoint a small value of the distribution parameter B is 
desirable. Again it can easily be seen that, consistent with the 
previous finding, it is desirable to use a low molecular weight gas. 
However, for given flight conditions, vehicle and injected gas, a low 
value of B implies a high injection mass-flow rate, and hence an 
increase in the weight of coolant required. Thus, it is apparent that 
the vehicle must be designed to keep the mass-flow rates low and still 
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keep B small. This can be attained to some extent through control of 
the shape parameters I s and R. Consider first the weighted surface-area 
integral I s . 

It is recalled that this integral is defined as follows: 


Is 




p . & 

%S R R 


Calculations have been made of the values of I s for hypersonic flows 
about hoth cylindrically r blunted wedges and spherically blunted cones. 
The results are shown in figures 5 an d 6, respectively. The theory of 
Lees (ref. 9) was used to obtain the heating- rate distribution for both 
planar and axisymmetric flows. The parameter b appearing in the 
planar case (i.e., fig. 5) is the span width. 

The results presented in figures 5 and. 6 tend to indicate that 
bodies with relatively small wedge or cone angles will have the lowest 
values of I s . In any direct comparison, however, it is necessary to 
hold some physical property of the bodies constant. The following 
discussion illustrates how such comparisons might be made. 

It has been shown (see, e.g., eq. (26b)) that the required mass of 
film coolant for various types of operation is proportional to: 



Now, consider for example an atmospheric entry at given entrance velocity, 
lift-drag ratio, and entrance angle. It can be shown that 


q A,s _ ( W Y /2 
AH A * \CjykRj 

or for a given over-all vehicle density 

^A^s f Vol 

AH A VciW 


(33) 


Lees has also shown in reference 10 that the distribution for planar 
flow differs from that for axisymmetric flow by only a few percent. 
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It is nov possible to define a drag-surface integral under hypersonic 
flight conditions so that 


Cd a s (5i>: 


(34) 


vhere 


Ip 1 ^ 5 sin6i S 

and where 5 is the local slope of the "body. Calculations have been 
made of the variation of Ip with s/R for spherically blunted cones. 
Pressure coefficients obtained with Newtonian Impact theory were used 
(i.e., Cp = 2 sin 2 6) and the results are shown in figure 7. 

From equations (32), (33), and. (3*0 It is possible to write the 
following : 


/, r -,\ 1/2-1/ 2 h 

m F * (Vol) R — - 

X P 

Thus for vehicles of the same volume, the one with the lowest value of 
R 1/ ' 2 (l g /lp 1 / 2 ) will require the lowest weight of coolant. 

From the information supplied by figures 6 and 7, the variation of 
the parameter R 1 / 2 (l s /lp 1/ ’ 2 ) with nose radius R has been calculated 
for spherically blunted cones with a volume equal to 5 0 cubic feet. The 
results are shown in figure 8, The curve labeled "Spherical Segment" is 
the limiting envelope where all the volume (50 cu ft) is contained in 
the spherical segment of the nose. For this example case at least, it 
is indicated that certain body shapes require less coolant than others. 
In particular, for bodies with small nose radii, small cone angles 
appear to have advantages for use with film cooling. 


Illustrative Examples 


To illustrate more graphically the effect of film cooling on the 
heat transfer to blunt bodies and to show additionally the relative 
importance of the parameters derived in the analysis, several example 
cases have been calculated. 
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F-tyecL flight conditions .- In the first example, a blunted cone with 
a half - an gle of 30 ° and a spherical nose of 1-foot radius was considered 
to he flying at an altitude of 200,000 feet and at a velocity of 16,000 
feet per second. The stagnation heating rate is then about 70 Btu/ft 2 -sec. 
The heating-rate distributions over the body for various amounts of 
helium film injection are shown in figure 9(a ) . In each case, the 
helium Is assumed to be injected at or ve ry near the stagnation point 
with an injection parameter, (l - -V Mp/MA. AHp 0 /AHa) = 1. It should be 
remembered that these curves are for very particular flight conditions 
and body shape and are only indicative of the distributions obtained 
with film cooling. With this restriction in mind, it is noted that 
within the framework of the analysis, the injection of sufficient helium 
can have a very pronounced effect on the heating- rate distribution and 
that the cooling may be effective several nose radii downstream of the 
point of injection. 

The effect on t he heat ing-rate distribution of changes in the injec- 
tion parameter (l - Mf /M ft AHf ; o/AHa) can "b e seen in figure 9(t>) • Here 

the same flight conditions and body shape are assumed as in the previous 
case shown in figure 9(a). In this case, however, a constant mass rate 
of injection of helium is assumed at 10” 3 slug/sec. Heating distributions 
with various values of the injection parameter are shown. It can be seen 
that the heating-rate distribu tion r eturns to that with no injection as 
the injection parameter (l - si Mp/M^ AHp 0 /aH^) approaches zero (i.e., as 
AHp^o is increased) . 

Nose radius is also an important parameter, as discussed earlier. 
Changes in nose radius produce two effects on the distribution parameter 
B. An increase in nose radius increases the surface area but also 
decreases the stagnation point heating rate. Some of these effects of 
nose radius on the heating-rate distribution are illustrated in figure 
9(c) . Again the same flight conditions were used as in the previous 
examples. An injection parameter of 1.0 and a mass rate of helium 
injection of 10“ 3 slugs/sec are also assumed. Results for various nose 
radii are shown. The injection has the largest effect for the bodies 
with the smallest nose radii. In this case, the primary effect is that 
(with the larger nose radii) a given mass rate of injection must cover 
a larger surface area and is thus less effective downstream of the 
injection point. 

Application to atmosphere entry .- Other examples have also been 
calculated to evaluate more clearly the use of film cooling for heat 
protection during atmosphere entry. The vehicle considered in these 
calculations is shaped similar to the Mercury capsule and is shown in 
figure 10. In these examples two entry trajectories for nonlifting 
vehicles, one from circular satellite orbit, and one from a parabolic 
orbit with a maximum deceleration of 8 g were considered. The vehicle 
characteristics used are summarized in the following table. 
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Entry 


W 


R W/CpA I s 


Circular 2500 lb 7 ft 

Parabolic 4000 lb 10 ft 


36.5 lb/ ft 2 

28.5 lb/ft 2 


0.77 

.77 


The normalized heating rate and velocity histories for the entries under 
consideration were obtained with the aid of reference 11 and are shown 
in figure 11, Only front face heating will be considered. 


It is suggested that film cooling should be used in connection with 
and as a complement to other means of heat protection. In this way the 
flexibility and possibly the efficiency of the over-all heat protection 
system can be improved. To illustrate this point consider a beryllium 
heat sink protecting the described vehicles entering along the described 
trajectories. The approximate total weights of heat-sink material 
required for protection under these conditions are shown in the following 
table: 


Entry Weight of heat sink; lb 

Circular 3^0 

Parabolic 1050 


The injection of helium can be used in conjunction with a heat 
sink of reduced weight. In the present examples It was considered that 
the helium was injected in such a way as to keep the integrated heating 
rate accepted by the vehicle at one-half the value existing without 
injection during the entire entry. Thus the total heat load for an 
entry is reduced by one-half. From figure 4, this reduction requires, 
for an injection parameter of unity, a distributional parameter, B, of 
1.60. Thus the mass of helium required is 


m^ 



*A,B 

2 % 


dt 


mp 




dt 


Since the entry heat load is reduced by one-half, the amount of 
heat -sink material required is also reduced by one-half. Accordingly, 
the results of the calculations for the example trajectories give the 
following total weights for the heat-protection system. 
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Entry 

Circular 

Parabolic 


Weight of 
heat sink 

180 lb 
525 lb 


Weight of 
helium 


12.5 lb 

22.5 lb 


The apparent veight savings due to helium injection are not as 
large as might be indicated by these numbers. Additional veight is 
required by the system for tankage and for flow regulation vhich may 
be as high as an order of magnitude greater than the helium veight. 


CONCLUDING REMARKS 


A simple method for the estimation of heating rates to blunt shapes 
vith film cooling has been developed. The method indicates that the 
primary requirement of film cooling is the replacement of the high- 
enthalpy boundary layer vith a lov-enthalpy film. As this film flovs 
along a body, the heat convected to the body increases somevhat as a 
result of an increase in enthalpy from the absorption of part of the 
convective heat transfer by the film. Results obtained vith the method 
vere compared vith existing experimental data. Unfortunately, no real 
conclusions could be made as to the accuracy of the analysis in predicting 
the cooling effect of film injection. The method vas also used to indicate 
some of the important parameters for the attractive performance of a 
film-cooling system. It vas found that (l) the molecular veight of the 
film gas should be lov, and (2) the enthalpy potential of the film 
should be small. 

Application of the analysis vas also made to the calculation of the 
veight of a heat protection system utilizing helium film cooling in 
combination vith a beryllium heat sink. This calculation vas made for a 
vehicle shaped after the Mercury capsule concept and for nonlifting 
entries from circular and parabolic speeds. The results of these calcu- 
lations vere compared to the veight s of a simple, all-beryllium heat 
sink. A significant veight saving vas indicated. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., March 28 , 1961 
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Figure 2.- Comparison of predictions with present theory with existing 
experimental results (nitrogen infection) . 




Data (ref. 6) 
Air injection 



Figure 3 •- Comparison of predictions with present theory with existing 
experimental results (air injection) . 
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Figure 5.- Weighted surface area integral for blunted wedges. 
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Figure 6.- Weighted 


surface area integral for “blunted cones. 
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Figure 7*- Drag surface integral for Hunted cones 
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(c) Effect of nose radii 


Figure 9 *“ Concluded. 
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